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O ■ ABSTRACT 

(N 

Context. Lyman Break Galaxies (LBGs) are widely thought to be prototypical young galaxies in the early universe, particularly 
representative of those undergoing massive events of star formation. Therefore, LBGs should produce significant amounts of X-ray 
emission. 

Aims. We aim to trace the X-ray luminosity of Lyman Break Galaxies across cosmic time and from that derive constraints on their 
star formation history. 

Methods. We utilize the newly released 4 Ms mosaic obtained with the Chandra X-ray Observatory, the deepest X-ray image to date, 
alongside with the superb spectroscopic data sets available in the CDF-S survey region to construct large but nearly uncontaminated 
samples of LBGs across a wide range of redshift (0.5 < z < 4.5) which can be used as input samples for stacking experiments. This 
approach allows us to trace the X-ray emission of Lyman Break Galaxies to even lower, previously unreachable, flux density limits 
Q^) ■ (~ 10~ 18 mWirT 2 ) and therefore to larger redshifts. 

Results. We reliably detect soft-band X-ray emission from all our input redshift bins except for the highest redshift (z ~ 4) one. 
From that we derive rest-frame 2-10keV luminosities and infer star formation rates and stellar masses. We find that star formation in 
LBGs peaks at a redshift of z pea k ~ 3.5 and then decreases quickly. We also see a characteristic peak in the specific star formation rate 
(sSFR=SFR/M„) at this redshift. Furthermore, we calculate the contribution of LBGs to the total cosmic star formation rate density 
j.^ ■ (SFRD) and find that the contribution of LBGs is negligible. Therefore, we conclude that most of the star formation in the early 

universe takes place in lower luminosity galaxies as suggested by hierarchical structure formation models. 
Conclusions. 
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OO ■ 1. Introduction I2002I) . Therefore, large samples of LBGs across the entire red 
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Lyman Break Galaxies (LBGs) are largely considered as the 2009; B asu-Zvch eUd1l201 ll) to z > 7 (iBouwens eTallfeoiO. 



shift r ange from z _ 1 (fBurgarella et alJl2007l:lHaberzett 



. bright end of the distribution of normal star-forming galaxies 201 ll |Stok^|201_iJ) are publicly available, typically com- 

■ across a wide range of redshift.They exhibit a significant scat- prising hundreds of objects m a compr ehensive form. Even at the 

I ter in terms of mass (Mannucci et al. 2009), with both a pro- highest reds hifts, considerable work is being done utilizing fore- 

(N . nounced low-mass ([Weatherley & Warren||2003|) and high-mass front e q U ip me nt such as the new Wide-Field Camera 3 (WFC 3) 

■ d BarmbyetalJ |200j) fraction among the entire LBG popula- aboard the HubMe Space Telescope (HST) . Therefore, LBG can 
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tion. Since they are - thanks to their exceptionally strong star didat es with redshifts as large as z = 10 are being discovered 
formation activity - relatively easy to select at various red- dBouwens et alJl2011^h . but because of their extreme faintness 



shift ranges from photometric observations only, they have been (Hab = 2g 9) lacking spectrosco pic confirmation which has only 

H . playing a key role in galaxy evolution studies all over the last been done for LBGs up to z ~ 7 (IVanzella et alJl20Tl . There is 

W ■ two decades. This is foremost due to the very effective "drop also cons i der a ble work chara cterizi ng LBGs in terms of their 

out" selection technique established in the earlyJ990s (see e.g. environment (Tasker & Brvan 2006; Cooke et al. 2010). 

ISteidel & Hamiltonlll993t ISteidel et al.lll996l Il999h . This tech- 

nique utilizes the strong absorption of all light emitted blue- Since LBGs are traditionally selected in the optical wave- 
ward the rest-frame wavelength of the Lyman limit at 912 A. length regime, recently extending to the near-infrared (NIR) 
Because this produces a very pronounced step in the typical as the WFC 3 aboard HST became operational, ancillary ob- 
spectrum of an LBG, one can select such objects by searching servations are necessary to characterize these objects over the 
deep imaging data for sources detected only in longer wave- entire electromagnetic spectrum. This has been done by sev- 
length filters but not in short wavelength filters. This behavior eral groups, successfully detecting individual LBGs at moder- 
of a source "dropping out" from being detected below a cer- ate redshifts (z < 3) also in the mid- and far-infrared (MIR 
tain wavelength is nowadays the major tool for selecting can- & FIR) regimes (iRigopoulou et alJ 120061 : iMagdis et al . 2008; 
didate high-redshift sources. Since its first application in the Burgarella et alJ 1201 ll) thanks to other space-based facilities 
early 1990s, the technique underwent a massive d evelopment such as the Spitzer Space Telescope and, recently, the Herschel 
and heavy usage troughout the community (review bv lGiavaliscol Space Observatory. At even longer wavelength, ground-based 
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observations in the sub-millimeter and radio regimes have been 
conducted, yieldin g only a few detectio ns of individual objects 
at 850/itn (e.g. Cha pman & Casevll2009f) whereas there is no di- 
rect detection of a significantly redshifted LBG to date except for 
strongly lensed sys tems such as the "8 o'clock arc', for instance 
(IVolinoetal.1120101) . 

To access even more redshifted LBGs over the entire spec- 
trum, stacking techniques have been successfully applied by 
many groups. For instance, the radio properties of LBGs have 
been discussed wi th respect to their star formation activity by 
ICarilli et ail d2008h . utilizing the deep Very Large Array (VLA) 
1.4 GHz (or 20 cm) observations of the COSMOS field. Similar 
stacking investigations of the star formation history of a sample 
of ver y high redshift (z > 7) LBGs were done by Labb e et al.l 
(2010) in the NIR regime, utilizing ultra-deep Spitzer data in 
the Chandra Deep Field South (CDF-S). This particular sur- 
vey field is also extremely valuable to study LBGs at the short- 
est wavelength since it comprises the deepest X-ray observa- 
tions obtained so far. However, with less deep data significant 
work on stacking LBGs in the X-ray regi me has already been 
done. A first attempt has been done already Bra ndt et al.1 {2001) 
in the Chandra Deep Field North (CDF-N) with a 1 Ms expo- 
sure. Stacking a sample of only 24 LBGs with redshifts be- 
tween 2 < z < 4, they found a soft-band signal at a signif- 
icance level of 99.9%. From that, they calculated an average 
X-ray luminosity of 3.2 10 41 erg s" 1 in the rest-frame 2-8 keV 
band, comparable to th e most luminous local starburst galax- 
ies such as NGC3256 (Mora n et al.lll999l) . A similar result was 
obtained by Nandr a et al.l (120021) who extended the small sam- 
ple of iBrandt et al.l to a statistically more robust number of 148 
LBGs. With this larger sample, they were able to exclude LBGS 
containing Active Galactic Nuclei (AGN) and hence giving an 
estimate for the star formation rate of an average LBG at z ~ 3 
of abo ut 60 M yr~ 1 . A next step was done by iLehmer et al.l 
(2005) utilizing Chandra data from both the Great Observatories 
Origins Deep Survey (GOODS) north and south fields with ex- 
posure times of 2 Ms and 1 Ms, respectively. With a large sam- 
ple of LBGs comprising a few thousand galaxies, they were 
able to do their analysis in different redshift bins in the range 
3 < z < 6, detecting a significant signal up to z = 4. They found 
the average star formation rate at this redshift decreases to 10- 
30 M yr" 1 , and even further to higher redshifts since they were 
not able to make statistically significant detections at any red- 
sh ift greater than z = 4. A complementary result was obtained 
bv lLairdet ai1 (f2006) who identified a sample of direct X-ray de- 
tections of LBGs at intermediate redshifts (z ~ 3), supporting the 
notion that LBGs are the "tip of the iceberg" in terms of star for- 
mation in their respective epoch. The statistical X-ray detection 
of high-red shift LBGs (z > 5) was attempted for the first time 
by ICowie etldl (1201 2|) utilizing photometrically selected sam- 
ples of LBGs comprising several thousand sources in the CDF- 
S. Although with such a large number of stacked sources, they 
were not able to detect X-rays from the high-z LBG population, 
placing an upper limit of 4 10 41 erg s" 1 on the X-ray luminosity 
in the rest-frame 4-15 keV band at z = 6.5. 

With this paper, we attempt to widen the redshift range of X- 
ray studies of LBGs down to z = 1 since space-based facilities 
such as GALEX and SWIFT now provide reliable LBG samples 
selected in the ultraviolet (UV). Furthermore, we want to trace 
the star formation activity and stellar mass build-up of the aver- 
age LBG from z = 1 onwards. 

Throughout this paper, we adopt a standard flat ACDM 
cosmology with Hp = 70kms~ 1 Mpc~ 1 and £2a = 0.73 
dKomatsu et al .11201 ll) . 



2. X-ray data and LBG selection 

In this paper, we use the 4 Ms Chandra Deep Field South data[3 
as basis for the X-ray stacking procedure. These data comprise 
mosaics of an area of 464.5 arcmin 2 in both the soft (0.5-2 keV) 
and the hard (2-8 keV) X-ray bands. With an effective exposure 
time of 3822 ks, this is the the deepest X-ray data set available 
to date, reaching an on-axis flux limit of 9.1 10~ 18 er g s" 1 cm" 2 
and 5 .5 10 _17 ergs 1 cm" 2 in the soft and hard band (iXue et al.l 
1201 lh . 

To create a parent list of Lyman Break Galaxies to be 
stacked, we use the extreme wealth of spectroscopic data avail- 
able in the CDF-S. To use only spectroscopically confirmed 
LBGs assures both a clean list with only marginal contamination 
by mis-classified objects, e.g. lower redshift interlopers, as well 
as accurate redshifts for the construction of several redshift bins. 
In particular, we use d the VLT/VIMOS spectr oscopic surveys 
dPopesso et al.l 120091) and Balestra et al. (2010) as wel l as the 
VLT/FORS survey b\ IVanzella et al.l (120051 120M |2008l) . These 
catalogs provide a detailed spectroscopic classification and are 
therefore ideal for reliably selecting LBGs. The main photomet- 
ric selection criterion was based on a (/- and B-band drop-out 
search utilizing a color cut U - B > 1.2 or B — V > 1.2, 
respectively, where the red part of the spectrum must stay flat 
(e.g. V - z < 1.2). The main spectral feature used in these sur- 
veys for deeming a source an LBG is the characteristic break 
in their spectrum blueward of the Lyman Limit at rest-frame 
912 A which was assessed by comparison to various template 
spectra of Lyman Break Galaxies and other, possibly contam- 
inating, source types. The spectral classifications provided by 
these surveys are very detailed, allowing e.g. for a split in Lya 
emitters and absorbers. However, in order not to bias our sam- 
ples in a certain way, we chose to not use this criterion but solely 
select LBGs (wether they are Ly a emitters or absorbers) in dif- 
ferent redshift bins. For a complete overview of our input list see 
Tab. Q] 

To extend our LBG sample to redshifts lower than z = 2.5, 
the approximate limit for the spectroscopic classification of 
LBGs using optical spectra, we had to take a closer look into 
the UV data of the CDF-S because at these low redshifts, the 
Lyman break occurs in the near UV range. Using the well- 
know n dropout technique to id entify Lyman Break galaxy candi- 
dates, iBasu^Z^ch^etaD |201l|) compiled a list of candidates with 
data from the SWIFT Gamma-ray observatory's UV/optical tele- 
scope (SWIFT-UVOT). The SWIFT satellite, although dedicated 
to high energy astrophysics, is equipped with a small (30 cm di- 
ameter) telescope sensitive to wavelength between 170/im and 
650 fim which observed an area of 266 arcmin 2 in the CDF-S for 
approximately 60 ks, reaching a limiting magnitude in the U- 
band of 24.5 AB-mag. These observations lead to a list of can- 
didate Lyman Break Galaxies with 43 objects, spanning the de- 
sired redshift range. This sample, dubbed BZ-UV, is also listed in 
Tab.Q] Because their selection mechanism is based only on pho- 
tometric data, we point out that the reliability of these objects is 
not comparable to the spectroscopic lists we use for higher red- 
shift objects. Therefore, and because there are much less LBG 
candidates identified via UV observations, we expect our stack- 
ing results in this redshift range to be significantly less robust 
than in the other redshift bins. 

We have taken special care in order to avoid any con- 
tamination by X-ray emission from a faint AGN that may be 



1 http : //www2 . astro . psu.edu/users/niel/ 
cdfs/cdfs-chandra.html 
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Table 1. LBG selection parent tables and redshift bins. 



ID 


parent list 


redshift range 


median z 


age" [Gyr] 


# selected objects 


# stacked objects 


ref-lowZ 


Wuvts et al. (2008) 


0.5 <z< 1.0 


0.68 


7.52 


143 


141 


BZ-UV 


Basu-Zvch et al. (201 D 


1.0<z<2.0 


1.33 


4.82 


43 


35 


Spec-Ul 


VIMOS & FORS 6 


2.5 < z < 3.0 


2.67 


2.49 


201 


201 


Spec-U2 


VIMOS & FORS 6 


3.0<z<3.5 


3.33 


1.95 


143 


92 


Spec-Bl 


VIMOS & FORS'' 


3.5 < z < 4.0 


3.70 


1.73 


63 


19 


Spec-B2 


VIMOS & FORS 6 


4.0 < z < 4.5 


4.16 


1.50 


39 


34 



Note s. (n) Time sinc e the Big Bang accordi ng to the cosmology defined in Sect.[JJevaluated at median redshift of the bin. (6) Popesso et al. (2009), 
iBalestra et ail i201d) . IVanzella et al] d2008h . 



hosted by our stacked sources. The rejection of possible AGN 
hosts was done in four stages: (i) The optical spectra of all 
sources in the input lists were examined both by eye and ac- 
cordin g to the classical BP T diagnostic diagram for AGN ac- 
tivity (Ba ldwin et alj IT98lT) . Any source with a marginal sign 
for AGN activity was discarded from the list, (ii) We used 
the deep Spitzer photometry in the CDF-S (GOODS project, 
see http : //irsa . ipac . caltech . edu/data/GOODS/) to ex- 
amine our input LBGs in the mid-infrared color -color diagram 
according to the criteria from IStern et al.l (2005) and excluded 
all sources lying in the AGN regime of this diagnostic dia- 
gram, (iii) We cross-correlated our input so urce lists to the Very 
Large Array (VLA) 20 cm observations bv lMiller et aLI (l2008h 
and plotted them against the Spitzer/MlPS 24 fim flux densities 
to obtain a measure for the radio-infrared correlation amongst 
our source s. Such sources t hat deviate from this correlation as 
defined in iMao et ail d201 ll) by showing a more than a 5 -fold 
excess in radio emission most likely originating from an AGN 
were discarded. 

In order to compare our findings for star formation rates 
which will be calculated from the stacked X-ray luminosities, 
we use data of our selected objects from other survey projects in 
the CDF-S. Primarily, the COMBO- 17 survey dWolf et al.ll2003l) 
is used for its broad range of data available over a large area. A 
deeper but more narrow view (in terms of covered a rea) is in- 
trodu ced by the use of data from the GEMS catalog dRix et alJ 
l2004t) with Hubble Space Tele scope (HST) obser vations in two 
filters and the FIREWORKS dWuvts et alJl2008l) catalog, con- 
taining deep data across a wide range of the electromagnetic 
spectrum (UV to MIR) but for only a fraction of our selected 
LBGs. In particular, we used the optical and NIR ancillary data 
to infer parameters such as stellar mass for all objects in our 
stacks. Furthermore, FIREWORKS was used to construct a low 
redshift control sample of galaxies. This sample should cover the 
very low redshift end in our analysis and be consistent of galax- 
ies of all kinds to obtain a completely unbiased view of X-ray 
emission from typical galaxies to be later compared to our LBG 
results. To compile this reference sample (dubbed ref-lowZ in 
Tab. [TJ), we only applied a redshift cut of 0.5 < z < 1 .0 to all 
objects classified as galaxy in the COMBO-17 catalog that also 
have FIREWORKS counterparts with reliable spectroscopic red- 
shifts. This left us with a fairly large sample of more than 200 
objects, therefore we expect the stacking to deliver an unambigu- 
ous X-ray detection. 



3. The stacking procedure 

Before starting the stacking process for the various source sam- 
ples defined above, special care was taken to exclude individu- 
ally detected sources from the samples. Therefore, in a first step, 



we cross-matched all our sources as summarized in Tab.Q]to the 
Chandra 4 Ms source catalog bv lXue et alj 12011 and excluded 
the 13 sources that have X-ray counterparts in this catalog which 
are mostly lower-z objects. For all remaining sources, 10" x 10" 
(resp. 20 x 20 pixels) cutout images from the 4 Ms soft-band mo- 
saic centered on the (optical or UV derived) source position were 
created and inspected by eye whether there is a clear but uncat- 
alogued X-ray counterpart. This was only true for two sources 
which were then also discarded from the further process. 

This left us with six source lists (according to our six red- 
shift bins from Tab. [TJ containing only objects without unam- 
biguously detected X-ray counterparts to start the actual stacking 
algorithm with. 

3.1. The stacking algorithm 

There has already been some eff ort in doing X-ray s tacking of 
faint sources as e.g. presented in iNandra et al.l {2002), and also 
stacking proce dures in other wa velength ranges, e.g. radio as 
described by Carilli et al. (2008), are known. A recent attempt 
al so using the new 4 M s Chandra CDF-S mosaic is presented 
in ICowie et alj d2012l) . They attempt to trace the X-ray emis- 
sion of galaxies out to z = 8 by utilizing high-z galaxy sam- 
ples compiled from new Hubble Space Telescope Wide Field 
Camera 3 (WFC 3) observations of the Hubble Ultra-Deep Field 
(HUDF) which is part of the CDF-S and applying a weighted 
mean stacking algorithm based on various quantities (off-axis 
angle, aperture radii for flux extraction, exact model of noise in 
aperture) not related to the sources themselves but to the techni- 
cal layout of the X-ray observations and the reduction process. 
This technique may enhance the resultant S/N ratio of the fi- 
nal stacked images, it bears the risk of introducing biases or 
unconstrained statistical effects just because of its complexity. 
This is discussed in m o re det ail e.g. in lLehmer et alj d2005l) and 
iHickox & Markevitchl (j2()07). Therefore, we chose to apply a 
much simpler and straightforward stacking method that, on the 
one hand, guarantees to not introduce any biases but for the even- 
tual cost of S/N in the final stacked images. Our stacking routine 
works as follows. 

We first consider the typical angular size of a galaxy at var- 
ious redshifts to determine a reasonable size for the images to 
be stacked. Since within our adopted cosmology the line ar scale 
at eve ry redshift z > 1 is nearly constant with 7kpc/" (Wright 
2006), the largest g alaxies (determined from their optical mor- 
phologies as e.g. bv lTrujillo et~ai]l2006l) at these redshifts with a 
diameter of ~ 10 kpc would stretch over about 1 .5" or three pix- 
els given the pixel scale of Chandra 's ACIS detector. Therefore, 
we chose to extract fluxes from a 5 pixel diameter aperture to 
also account for astrometric inconsistencies between the parent 
galaxy catalogs and the 4 Ms mosaic. Note that this astrophysical 
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Fig. 1. The local Lyman Break Galaxy analog VV 1 14 as seen by 
Chandra (color co mposite of soft- and hard-band images from 
Grim es et alj|20 06). The overlaid grid indicates the pixel size of 
the ACIS detector if VV 1 14 would be at z = 3 and our 20 x 20 
pixels stack-image layout. 



motivation of an aperture siz e is very much diffe r ent from other 
stacki ng studies such as e.g. lNandra et alj d2002); Le hmer et alj 
(120051) . They defined their aperture sizes empirically by testing 
various diameters and then adopting the one yielding the best 
S/N ratio. Despite the fact that those two methods are fundamen- 
tally different, the final aperture diamters both our astrophysical 
and the empirical approach yield are in very good agreement. 
To finally get a measure of the background around every source, 
we chose the stacking images to be 20 by 20 pixels, four times 
the aperture diameter. The layout of our images that go into the 
actual stacking procedure is shown in Fig.Q] 

Before actually stacking these 20 x 20 pixels images, our 
algorithm sorts out all images that (i) show an integrated flux 
within the pre-defined 5 pixel aperture (centered on the image 
center) that exceeds three times the median absolute deviation 
(MAD) of the noise in the remaining part of the image and (ii) 
all images with a total flux across the entire image that is zero. 
The first step in our algorithm is to reject all marginally detected 
sources to avoid any contamin ation of su c h sour ces that are nei- 
ther in the source catalog by IXue et a D (1201 lh nor visible by 
eye. This can be due to e.g. AGN which are faint at optical and 
infrared wavelength and t herefore are not visible the FORS or 
VIMOS spectra (see e.g. lZinn etaDl201 lh . Since at a redshift 
Z — 3 the X-ray luminosity necessary to produce an actual detec- 
tion in the 4 Ms mosaic is equivalent to a star formation rate of 
about 200 M Q yr _1 , an actual X-ray detection is a strong hint for 
additional AGN activity because LBGs are in general not har- 
boring such extreme starbursts. The second criterion is applied 
because the 4 Ms mosaic is smaller than the region for which 
spectroscopic observations have been carried out. Hence, there 
may be sources in our list that do not have X-ray coverage which 
should of course be sorted out in order not to contaminate the fi- 
nal stack. 

To also exclude X-ray detected sources that are just by 
chance within our stacking region, we define a second aperture 
of 10 pixels diameter again centered on the image center. The 
remaining part of the image is used to determine the background 
noise level. If the total flux measured within this 10 pixel aper- 
ture indicates a detection with a significance of more than 90% 
with respect to the measured local background, it is rejected for 
the actual stacking process. This second rejection step is neces- 
sary to remove X-ray sources that are close to the stacking loca- 



tion but are not associated with the LBG that should go into the 
stack. Because those unrelated sources could add flux to the fi- 
nal stack which is not physically related to the LBGs we want to 
examine, this second rejection step is crucial in terms of a clean 
stacking sample. Tab. Q] summarizes the number of sources se- 
lected for the six input coordinate lists and the number of sources 
actually stacked after this rejection step. 

The actual stacking then consists of a simple average stack- 
ing without any weighting in order to not introduce potential 
biases due to more complex statistical treatments. We point out 
that, as common in most other wavelength ranges where the un- 
derlying noise distribution is Gaussian, we cannot use a median 
stacking because most of the pixel values stacked are zero, hence 
the median would always be zero. This is particularly unfor- 
tunate because the median is known to be much more robust 
against outliers (e.g. one unusually bright pixel that immediately 
increases the average would just be ignored by the median). This 
is the main reason for applying the rejection as outlined above. 

The flux in the six final stacked images is then extracted 
by integrating over the pre-defined 5 pixel aperture. To max- 
imize the flux covered by the aperture, it is actively centered 
on the brightest pixel in the image and not fixed at the image 
center. To correct for noise within the aperture, the background 
noise is calculated from the remaining part of the image (every 
pixel not covered by the flux extraction aperture) and then sub- 
tracted from the raw flux measured within the aperture. Because 
not every pixel in the 4 Ms mosaic has exactly the same expose 
time, we used the exposure maps distributed with the actual sci- 
ence frames to determine a precise total exposure time for every 
stack. By dividing the noise-corrected number of counts within 
our aperture by this total exposure time, we get the final X-ray 
flux in the respective band in counts per second. The conver- 
sion to physical units (here ergs" 1 cnr 2 =mWrrT 2 ) is then done 
by extrapolating a conversion factor between counts per secon d 
and flux based on the 4 Ms source catalog bv lXue et al.l d201 lh . 
To quantify the significance of detection in our final stacked im- 
ages, the total number of counts within our extraction aperture, 
dubbed x, is compared to the noise level of the stacked image, 
dubbed A. Assuming a Poissonian noise distribution, this noise 
level should be the expactation value of the distribution. Hence, 
we calculate a cumulated probability P which is the probability 
of measuring x or less counts where we expect A counts. If this 
probability exceeds 95%, we assume the stacked detection to be 
real. We point out that, in most cases, our detections are well 
above the 99% confidence level. To give an error on our mea- 
sured flux densities, we compute the number of counts necessary 
for a 68% (corresponding to lcr confidence level in a Gaussian 
distribution) signal and attribute this quantity to be the error of 
our measured flux. 

3.2. Testing the algorithm 

Considering the (mathematical) simplicity of our algorithm, we 
expect to avoid any biases due to statistical effects. Nevertheless, 
we ran the following tests on artificial data to ensure the perfor- 
mance of our stacking procedure in terms of reliability of stacked 
detections. Note that the main assumption of stacking - the noise 
goes down with the square-root of the number of stacked objects 
- is also valid for Poissonian noise distributions since they also 
exhibit (on the scales of the stacked sources) uncorrelated noise: 
Therefore, we artificially created images with 500 by 500 
pixels and a Poissonian noise distribution (as it is the case for X- 
ray images). Every image contains 27 point sources according to 
the PSF of Chandra at random off-axis angels and signal to noise 
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Table 2. Summary of tests of the stacking algorithm. The values 
given in the table are the median S/N ratios of the final stacked 
sources over all 100 test runs. 







stacking offset 




none 


5 pixels random 


0.1 


nd 


nd nd 


initial S/N 1.0 


2.91 


2.88 nd 


3.0 


8.88 


8.89 nd 



Notes, "nd" indicates that in the particular stack no source is detected 
(in the case for a stack of S/N=3 sources no source with a significantly 
higher S/N). 



ratios of 3, 1 and 0.1 to obtain nine barely detectable sources, 
nine not detected sources and nine sources well beyond the de- 
tection limit. Because for the Chandra 4 Ms mosaic, spec ial care 
was laid on astrometric calibration (see IXue et alj|201 ll section 
4.2) resulting in a median posit i onal o ffset between X-ray and 
radio sources from Mill er et aD d2008l) of 0.24", so clearly be- 
low the pixel size of Chandra which is 0.5", we ignored effects 
caused positional uncertainties and worked with the exact posi- 
tions of the artificial sources. 

We then stacked the sources according to their attributed S/N 
ratio using our algorithm as previously described. As usually as- 
sumed in stacking procedures, the noise level should decrease 
with the square root of the number of stacked objects. Therefore 
we expect the original S/N ratios to be increased by factor of 
three, making the barely detected ones well detected, the not de- 
tected ones barely detected and the sources well beyond detec- 
tion should remain undetected. These expectations were met in 
all three cases: the stack of the nine S/N=3 sources contains one 
nicely centered source for which we measured a new S/N ratio 
of 8.8, the nine S/N=l sources were stacked into a single source 
with S/N=2.9 and the nine S/N =0.1 sources were not visible in 
the stacked image at all. The fact that the S/N ratios reached in 
this test are always a bit below the theoretical expectation of the 
V^V-law is mainly just due to small number statistics. 

To further verify the algorithm and exclude random effects, 
e.g. from the fact that there may be a pixel with a high pixel 
value due to the long tail of the Poissonian noise distribution 
right at the site of an artificial source and hence unexpectedly 
increase their S/N ratio, we injected the same 21 sources to the 
same image again but with a constant offset of five pixels. The 
sources were then stacked again, resulting into stacked images 
nearly identical in S/N ratios to the first ones without offsets. 
For the final test, we stacked the sources with random offsets 
of up to 5 pixels. As expected, no source could be detected in 
these stacks (except of course for the nine stacked sources with 
an initial S/N=3). 

These three tests were repeated on 100 images to assure 
statistical significance. A summary of the tests is presented in 
Tab. [2] giving the median S/N ratios of the sources in the final 
stacked images (in case there is a source detected). 

To get a feeling of the stacking as an entire procedure, we 
also show a panel illustrating our routine (Fig. [2]): The left im- 
age shows an artificial source with S/N=5 sitting in a Poissonian 
noise distribution. The image in the middle is of the same size 
and noise level but with a source of S/N=l in the image cen- 
ter and the right image shows the resultant stack of nine such 
S/N=l sources (smoothed with a Gaussian of three pixel full- 
width-half-maximum) using our stacking algorithm. We mea- 
sured S/N=2.95 for this stack within a circular aperture of 5 pix- 



els diameter. Compare to Fig. [3] showing our stacks with real 
data. 

4. Stacking results 

Our final stacked soft-band images for all six input lists (resp. 
redshift bins) are shown in Fig. [3] a table summarizing the ex- 
tracted fluxes and detection probabilities is presented in Table[3] 

The flux extraction process was performed in two steps: 
First, all stacks were visually inspected for whether there is a de- 
tection or not. To quantify this, we then extracted the number of 
counts both within our pre-defined 5 pixel aperture (where the 
aperture was centered on the pixel with the highest count rate 
for the images where the visual inspection suggested a detec- 
tion, for the others the aperture was centered on the image cen- 
ter) as well as for the remaining part of the stacked image. From 
the distribution of counts outside the aperture, we constructed 
the noise of the stacked image by fitting a Poissonian distribu- 
tion to the histogram. This was done to verify the Poissonian 
behavior of the noise in general and to obtain the expectation 
value, A, for every image. From this expectation value, we es- 
timated the noise count rate that would be comprised by our 5 
pixel aperture and subtracted it from the actual count rate mea- 
sured within this aperture (correcting for pixels of which only 
a fraction lies within the aperture), yielding a count rate of the 
stacked sources only. Using our fit to the noise distribution, we 
then calculated the probability that such a count rate occurs just 
by chance in the respective image. Its inverse probability (so 1 - 
"by-chance-probability") is then attributed the confidence of de- 
tection. Adopting a common confidence limit of 95%, we con- 
sider all stacks with a confidence greater than 0.95 as detections. 
The source count rates (the ones cleaned from noise counts) and 
the confidence values are given in Table [3] 

We then converted the source count rates to physical flux 
densities assuming a power-law X-ray spectrum with a photon 
index typ ical for star-formi ng galaxies of F=2.0 as empirically 
found by iPtaket al.l d 1999b. see also the discussion on X-ray 
spectra of LBGs in lLaird et ail (H006). We also corrected for 
Galacti c extinction in the C DF-S (N m = 6.8 10 19 cirT 2 ) fol- 
lowing iKalberla et all J2005I) . To test our flux calibration, we 
compa red it to the original 4 Ms source catalog by IXue et ail 
(2011) finding them in very good agreement with their fluxes. 
Errors were assigned to all measured fluxes by taking into ac- 
count the image noise (more precisely, the 68% uncertainty of 
the extracted flux, dubbed 5f) and a possible 10% error due to 
the aperture size (dubbed <x). This error accounts for the fact 
that we chose our aperture to have a fixed size in contrast to 
other work where the aperture is empirically determined t o max - 
imize the S / N rati o for each stack, see e.g. lLehmer et alj d2005l) : 
lLaird et al.ld2006h : 

A/ = ^(6f) 2 + (cr-f)\ (1) 

where / is the measured flux of the detected source. 

To further validate our results regarding the X-ray luminos- 
ity computed assuming a power-law with photon index F=2.0, 
we calculated the hardness ratios (or upper limits, respectively) 
for each stack and check whether they are consistent with X-ray 
emission following such a power-law distribution. The hardness 
ratio is accordingly defined as HR — (H - S)/{H + S), where H 
and S are the counts in the hard- and soft-band, respectively. 
Since obscured AGN at the considered redshifts would show 
rather flat X-ray spect ra with pho t on ind ices T ~ 1 .0 - 1 .4, corre- 
sponding to HR > O.O lPark et al.l (T2008). the computed hardness 
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Fig. 2. Plot illustrating the different steps of our stacking procedure: (a) shows an artificially created source of S/N=5 at the image 
center within a Poissonian noise distribution, (b) is the same as (a) but with a S/N=l source and (c) shows the final stacked image 
of nine such S/N=l sources (smoothed with a Gaussian of three pixel full-width-half-maximum) using our stacking algorithm. 



Table 3. LBG stacking results. 



ID 


median z 


soft 


conf 


hard 


conf 1 




HR 


soft flux 


soft flux err* 


hard flux 


hard flux err 6 






counts 




counts 








10-' 8 mWnr 2 


10- 18 mWnr 2 


10- 18 mWnr 2 


10- |8 mWm- 2 


ref-lowZ 


0.68 


153.83 


0.999 


60.01 


0.907 


< 


-0.44 


2.15 


0.37 


4.41 


u.l. 


BZ-UV 


1.33 


36.36 


0.993 


25.12 


0.850 


< 


-0.18 


2.05 


0.24 


7.44 


u.l. 


Spec-Ul 


2.67 


112.26 


0.999 


100.74 


0.971 




-0.06 


1.10 


0.23 


5.20 


1.88 


Spec-U2 


3.33 


39.91 


0.958 


14.54 


0.645 


< 


-0.47 


0.85 


0.19 


1.64 


u.l. 


Spec-Bl 


3.70 


21.52 


0.974 


16.83 


0.722 


< 


-0.12 


2.23 


0.69 


4.19 


u.l. 


Spec-B2 


4.16 


24.14 


0.660 


7.09 


0.626 






1.40 


u.l. 


2.16 


u.l. 



Notes. For the convenience of the reader, we note that 1 mWm 2 =1 ergs 1 cm 2 . (n) Probability that this detection is real, assuming a Poissonian 
noise distribution. See text for further explanation. (A) u.l. indicates upper limit since no actual detection could be made (confidence < 0.95). 



ratios from Table[3]reveal that also the final stacks are completely 
consistent with X-ray emission originating from star formation 
rather than from AGN activity. For completeness, we note that 
also the soft-band detections and hard-band non-detections for 
the two low-redshift control samples are consistent with the X- 
ray luminosity derived from the soft-band fluxes and a power- 
law SED with r = 2. 

5. Astrophysical interpretation 

Table|4]summarizes our derived parameters that are discussed in 
the next sections. 

To obtain the rest-frame 2-10keV luminosity, we used our 
soft-band stacked fluxes and the corresponding median redshift 
of each bin. The median redshifts were calculated from the orig- 
inal spectroscopic redshifts as given in th e inp ut catalogs by 
Popesso et all (120091) . iBalestra et all (l2010h . and IVanzella et all 
2008) for our spectroscopic s amples or, respect i vely, f rom the 
photometric redshifts given in iBasu- Zvch et al.1 d201ll) for the 
BZ-UV sample or from the spectroscopic redshifts provided by 
IWuvts et"al] (12008) for the low-z reference sample comprised 
from the FIREWORKS catalog. To obtain a 2-10keV luminos- 
ity, we used the Portable, Interactive Multi-Mission Simulator 
(PIMMS) version 4.40 to get a 2-10keV flux based on an ex- 
trapolation of our stacked soft-band fluxes and a photon inde x 
T=2.0 typical for LBGs as e.g. found by lLaird et al.1 (120061) , 
This factor is adopted throughout the stacking literature, see 

2 http://asc.harvard.edu/toolkit/pimms.jsp 



e.g. iNandra et all d2002l) or lLehmer etal.1 d2005l) . Note that the 
stacked fluxes from Table [3] are already c orrected for Galactic 
extinction following Kalberl a et al.l ((2005). These fluxes could 
then be directly converted to luminosities with respect to the 
above specified cosmology. We did not account for any errors 
due to redshift uncertainties because the flux errors are by far 
dominating the luminosity errors. 

Star formation rates were calculated using the calibration by 
iRanalli et al.l d2003l) : 

SFR [Moyr- 1 ] = 2.0 x 10- 33 L 2 -iokeV (2) 

where L2-iokeV is in units of Watts. The corresponding errors 
take into account both the uncertainties i n luminosity a nd the 
intrinsic scatter of the calibration. Since Rana lli et al.l (2003) 
did not quantify this scatter, we computed it using their data, 
finding an rms scatter of roughly 20%. In addition to the in- 
trinsic scatter in that relation we point out that it is calibrated 
using a local galax y sample only and that recent studies (e.g. 
Diiks tra et al.ll2012b consider the possibility that the conversion 
factor between S FR and i2-iokeV may increase with redshift, 
introducing a potential underesti mation of star formation rates. 
The same conclusion is drawn by Svmeonidis et aT] (1201 lb who 
quantify their deviation to the IRanalli et al.l d2003l) calibration 
to be a factor of five. This w ould directly imply that, adopting 
the ISvmeonidi s et al.l d201 lb conversion factor, our star forma- 
tion rates would increase by a fact or of five. However , for further 
analysis we chose to stick to the IRanalli et al.l (120031) calibration 
just because it is the most widely used calibration throughout the 
literature. 
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Fig. 3. The final stacked soft-band images (all smoothed with a Gaussian of FWHM = 3 pixels for illustration purposes) used for 
further analysis: (a) ref-lowZ (P = 0.999), (b) BZ-UV (P = 0.993), (c) Spec-Ul (P = 0.999), (d) Spec-U2 (P = 0.958), (e) Spec-Bl 
(P = 0.974), (f) Spec-B2 (P = 0.560). The white circles indicate the 5 pixels diameter apertures used for flux extraction. 

Table 4. Derived parameters for our LBG samples. 



ID 


median z 


rest 2-10keV 


rest 2-10 keV err* 


SFR 


SFR err* 


M» 


M, err 


sSFR 


sSFR err* 






log(W) 


log(W) 


M Q yr" 1 


M Q yr' 


10'° M Q 


10 10 M Q 


Gyr" 1 


Gyr 1 


ref-lowZ 


0.68 


33.06 


32.30 


2.3 


0.40 


6.73 


1.61 


0.03 


0.01 


BZ-UV 


1.33 


33.65 


32.72 


8.93 


1.05 


3.35 


0.67 


0.27 


0.08 


Spec-Ul 


2.67 


33.97 


33.29 


18.67 


3.90 


1.16 


0.29 


1.61 


0.15 


Spec-U2 


3.33 


34.03 


33.38 


21.43 


4.79 


1.09 


0.28 


1.96 


0.80 


Spec-B 1 


3.70 


34.53 


34.02 


67.77 


20.97 


1.09 


0.21 


6.25 


1.49 


Spec-B2 


4.16 


34.42 


u.l. 


52.61 


u.l. 


1.15 


0.31 


4.56 


u.l. 



Notes. For the co nvenience of the reader, we note that 1 W= 10 7 erg s 1 . (n) Assuming a power-law X-ray spectrum wi th a photon index r=2. 
typical for LBGs dLaird et al.l2006h . This factor is adopted throughout the stacking literature, see e.g. lNandra et al.l d2002h or Lehmer et al. ( 2005). 



u.l. indicates upper limit since no actual detection could be made (confidence < 0.95). 



Since our final goal is to also quantify the specific star for- 
mation rate sSFR, defined as the star formation rate divided by 
the stellar mass of the galaxy, we need to get a handle on the 
stellar mass. To do so, we computed the mean /^-band magni- 
tude of the LBGs i n every stack using data from FIREWORKS 
dWuvts et al.l20 08). This mean A'-band magnitude was then con- 
verted to a stellar mass following the calibration by Dadd feTaT] 
(120041) : 



log(M t /10 n M o ) = -0.4(^-^10 



(3) 



where K is the total /f-band magnitude of the respective galaxy 
and K\\ = 20.15 (Vega system) the typical ^-magnitude of a 
10 M galaxy. A plot illustr ating the results for our LBG sam- 
ple is shown in Fig. [4] Since iDaddi et all (|2004) found this re- 
lationship for galaxies in the redshift range 1 < z < 3, we note 



that our low-z control sample may be affected by a systematic 
effects caused by its lower redshi f t. The scatter of this relation- 
ship is quantified by IDaddi et alj d2004l) to be ~0.1 dex, there- 
fore the error on M, takes into account both this scatter as well 
as the lcr uncertainty of the mean /T-band magnitude. The errors 
of the specific star formation rates are then calculated by sim- 
ple Gaussian error propagation. Note that, since the derived stel- 
lar masses for our LBG sample are just above the completeness 
limit, the higher-redshift bins might be affected by incomplete- 
ness issues leading to an overestimation of stellar mass. 



5. 1 . The X-ray luminosity of Lyman Break Galaxies 

Fig. [5] shows the measured X-ray fluxes versus redshift. In gen- 
eral, our results agree well with other X-ray stacking studies 
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Fig. 4. The stellar masses of our LBG sample with redshift. 
For comparison, the stellar masses derived from the entire 
FIREWORKS sample with spectroscopic redshifts are shown, 
too. The plotted completeness limit corresponds to a Ssigam lim- 
iting magnitude in the Zf s -band of 22.85 Vega-mag. 



such as, for example, Lehm er et alJ (|2005) who found a rest- 
frame 2-8 keV luminosity for a sample of 449 f/-band dropouts 
(z ~ 3) of (1.5 ± 0.3) 10 34 W and (1.4 ± 0.6) 10 34 W for a sam- 
ple of 395 bright B dropouts (z ~ 4), also seeing a nearly con- 
stant luminosity with increasing red shift. A similar result was 
also obtained bv lNandra et al.1 d2002h who give a 2-10keV rest- 
frame luminosity of (3.4 + 0.7) 10 34 W for their sample of 144 
[/-band dropouts in the Hubble Deep Field North. They also 
stacked a sample of 95 "Balmer Break" galaxies which are lo- 
cated at z ~ 1, comparable to our BZ-UV s ample with a median 
redshift of 1.33. For their z ~ 1 sample, INandra et all d2002l) 
report a 2-10keV rest-frame luminosity of (0.33 + 0.05) 10 34 W 
which also agrees w ell with our value from Tab. [4] Furthermore, 
INandra et al.l d2002) only find significant detections in the soft- 
band stacks just as in this work. The only hard-band stacking 
detection is made for the Spec-U2 bin (presumably because it is 
the bin containing the most stacked sources) with a hard-band 
flux of (5.20 + 1.88) 10 18 mWirr 2 . Together with the corre- 
sponding soft-band flux (see Tab. [3]), this gives a photon index 
T w 1.1, significantly lower than the value of 2.0 typical for star- 
formation activity. Since such low photon indices are more typi- 
cal for AGN, we assume that this one hard X-ray detection may 
be due to a small remaining contamination of an AGN which still 
survived the many steps of AGN rejection previously applied to 
all stacked sources. An alternative explanation may be the con- 
tamina tion by a strong e mission line at 7.47 keV as discussed 
e.g. in iFiore et all (120121) . However, since the effective area of 
Chandra becomes very small for such high energies, the contri- 
bution of photons with energies exceeding 7 keV only makes up 
a small fraction of the total counts in the 2-8 keV band. 

Ho wever, we highlight that recent studies by iLehmer et alJ 
(120121) hint a steep rise of number counts of normal star-forming 
galaxies just below the soft-band flux limit of the Chandra 
4 Ms mosaic. They report that at flux levels of the order of 
10~ I7 mWrrT 2 normal SFGs make up nearly half the X-ray 
number counts and that at even lower flux levels those galax- 
ies will completely dominate the population. Therefore, since 
our stacking analysis probes the faint X-ray population down to 



• stacked soft band 

♦ stacked hard band 

M 82 soft band 

4 Ms soft-band limit 




Fig. 5. The stacked fluxes of all samples in the observed-frame 
soft- and hard-band vs. redshift. Horizontal error bars indicate 
the width of the redshift bins, the solid line indicates the ob- 
served soft-band flu x of the low-z LBG analog VV 1 14 from 
iGrimes et all (120061) . VV 1 1 4 has an IR-derive d star formation 
rate of SFR = 48M yr _1 (ISoiferet al.lll989l) . The horizontal 
dashed line indicates the soft-band o n-axis detection li mit of the 
Chandra 4 Ms mosaic as reported bv lXue et alJ (1201 ll) . 



10 18 mWm 2 , our findings are also consistent with extrapola- 
tions from individually detected sources. 

5.2. Star formation in Lyman Break Galaxies 

From the rest 2-10keV luminosities of our various redshift bins, 
we estimated the star formation rat es for an average L yman 
Break Galaxy using the method by iRanalli et alJ 0003), The 
resulting SFRs are shown in Fig. [6] (left panel) alongside with 
several other SFR measurements for LBGs calculated from data 
taken at different wavelength ranges. Our X-ray derived SFRs 
peak at a redshift of about z pea k = 3.5 with the last measurement 
at z ~ 4 being an upper limit slightly below the peak value. 

Compared to SFR estimated from other wavelength ranges, 
we find at least a good agreement in the trend (increase to 
Zpeak - 3.5, then decrease)whilst the absolute numbers differ by 
a factor 2 to 5, depend ing on the particula r data set for com- 
parison. Especially the Karim et aT| d201 ll) values show much 
larger values which could be easily explained by the very dif- 
ferent selection criteria they used for their stacking input sam- 
ples: While we focus on LBGs, they selected their sample based 
on 3.6pm flux density to be able to split not only in redshift 
bins but also in stellar mass bins. However, the trend of increas- 
ing SFR to at least z = 3 is seen in both our LBG as well 
as their 3.6pm samples. To comp are with oth e r LBG samples, 
we plotted t he values obtained by ICarilli et alJ (120081) for lower 
redshifts and Smit et al ] (120121) for higher redshifts. ICarilli et al.l 
(2008) used LBGs in the COSMOS field to create U-, B- and V- 
dropout samples as input for a 1.4 GHz stacking analysis. They 
detected only the f/-dropouts (z ~ 3) in the radio stacks and 
find SFR = 31 + 1 M^yf 1 (adopting the SFR-L lAG nz cali- 
bration bv lYunet al.ll200ll) while the average UV-derived SFR 
for this sample is 17 M Q yr . These two values agree well with 
our X-ray based estimate of 21.43 + 4.79M yr~ 1 at this red- 
shift. The decrease with redshift is confirmed by the results of 
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Fig. 6. The star formation r ates (left) and the specific star formation rates (right) for differ ent LBG samples vs. r edshift. We compare 
our findings to the ones bv lKarim et al.l (1201 II) . iMagdis et al.l d20Tob . ICarilli et all d2008l) and lSmit etal] (1201 2|) using various other 
star formation rate indicators such as radio or UV luminosities. 



LSmit etal.l(l2012h who investigated the star formation in higher- 
redshift LBGs based on their rest-frame UV continuum slopes . 
Accounting for dust extinction following Meur er et al.l (jj999), 
they find SFRs of 35 M Q yr 1 and 24 M Q yr" 1 at redshifts of 3.8 
and 5.0, respectively. Comparing their results to other data, they 
see a peak in star formation rate between 3.0 < z pea k < 3.5 and 
a peak value of 50-60 M yr _1 , very comparable to our X-ray 
based findings. 

All together, we can confirm a peak in the star formation 
activity of LBGs at a redshift around z pe ak = 3.5 which is 
seen by various authors using various SFR estimation methods. 
However, the absolute values for star formation rate differ by 
quite a large margin which is either due to the different selection 
criteria for the galaxy samples investigated or inconsistencies in 
the calibration of the different SF R estimation methods in dif- 
ferent wavelength ranges, see e.g. iKurczvnski et al.l (1201 Oh or . 
This again shows that a uniform and consistent cross-calibration 
of SFR indicators is desperately needed, particularly in advance 
of upcoming large survey projects observing across the entire 
electromagnetic spectrum. 

Regarding the specific star formation rates (sSFR), a similar 
peaking trend with redshift is seen (Fig. [6] right pan el). Although 
our ga laxies are in general less massive than in the iKarim etaLl 
d201 ll) selection, the two samples are in better agreement while 
looking at the sSFR since it normalizes for the mass. Therefore 
it looks like the LBG sample is smoothly continuing the 3.6 ^m 
selection at higher redshifts. This suggests that at redshifts of 
z ~ 4, LBGs are typical for the galaxy population whereas at 
lower redshifts more massive galaxies are abundant. Since the 
upper limit obtained for z - 4.16 is significantly lower than 
the actual measurement at z — 3.70, a clear decrease for higher 
redshifts is also visible in sSFR. This peaking behavior of the 
specific star formation rate supports the widely adopted picture 
of stellar mass growth having a peak somewhere between red- 
shift 2 and 4, but contradic ts findings b y other authors who ob- 
served dFeulner et al.ll2005l) or modeled dKhochfar & Silkll201ll) 
the sSFR for a wide range of redshifts and find a more or less 
constant sSFR (at least for lower mass galaxies) from z — 2 
onwards to higher redshifts. An explanation for this difference 
could again be the di fferent sample selecti on criteria. While we 
looked at LBGs only, iFeulner et al.l (120051) did a more complex 



near-IR selection to again split in stellar mass bins. A compa- 
rab le LBG sample in th e Subaru Deep Field was investigated 
by lYoshida et all (2006) who also find a peaking sSFR with 
a peak value of about 0.1 Gyr -1 at z ~ 4. This is more than 
an order of magnitude lower than our peak value, but since 
they focus on more massive LBGs this difference is not a sur- 
prise. Th ey argue (according to the an alytical galaxy evolution 
model by iHernquist & S pringel 2003) that the trend of a peak- 
ing SFR/sSFR at redshifts between 3 and 4 can be explained 
by different parameters dominating the star formation process 
at various redshifts: At lower redshifts, the star formation ac- 
tivity is mostly governed by the cooling rate of the (molecular) 
gas residing in a dark matter halo whereas at higher redshifts 
the conversion of cold gas into stars is the dominating parame- 
ter regulating star formation. The transition between these two 
modes is marked by the peak of star formati on activity. This hy- 
pothesis is recently supported by the work of Re ddv et al.l (1201 2|) 
who also see a (mild) peak in the specific star formation rate of 
several Gyr -1 at z ~ 3. 

5.3. The contribution of LBGs to the cosmic star formation 
rate density 

Fig. Q shows the co-moving star formation rate density (SFRD) 
as derived from our star formation rates (Tab. [4]i and a cos- 
mology according to Sect. Q] As one can easily see, the peak- 
ing trend also continues in this plot, although the peak is much 
less pronounced. W e compare our estimates to values from 
lYoshida et al.l d2006l) and references therein (green diamonds in 
Fig. [7]) to find a similar peaking trend albeit at absolute SFRD 
values being about one order of magnitude higher than ours. 
These higher values can be explained by the d ifference in LBG 
sample selection between Yoshida et alJ d2006l) and our work (all 
our LBGs are spectroscopically confirmed whereas they employ 
a mixed spectroscopic/photometric selection process). 

However, despite the (relatively) small difference between 
the two LBG samples, there is a much larger difference when 
comparing SFRDs derived from LBG samples to estimates of 
the total cosmic star formation rate density as e.g. compiled by 
iHopkinsI d2004l) . The difference between total and LBG SFRD 
is about four orders of magnitude, hence deeming LBGs to be 
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Fig. 7. The contribution of LBGs to the total cosmic star for- 
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entirely negligible when investigating the bulk of star formation 
activity in the universe's history, even when considering poten- 
tial incompleteness of our LBG sample. This finding agrees very 
well with pr evious investigations of the cosmic star formation 
histor y (e.g. ISawicki &" Thompsonl l2006t iBouwens et al.ll20Q7l 
I2009h . From constraints on the UV luminosity functions at vari- 
ous redshifts, these authors drew conclusions for the question on 
whether the main fraction of star formation activity takes place 
at the luminous or faint end of the high-redshift galaxy popu- 
lation. Finding extremely steep faint-end slopes at all redshifts 
z > 2, they argue that most of the UV emission (and hence star 
formation) takes place in low-luminosity galaxies with only a 
few per cent contribution of the bright population. More quanti- 
tatively, ISawicki & Thompson! d2006l) derived characteristic lu- 
minosities (and the correspondin g absolute magnitude Al*) 
for LBGs by fitting ISchechterl (Tl 9761) luminosity functions (LFs) 
at redshifts 2.2, 3.0 and 4.0. They find a nearly constant value of 
Ai t = -21.0 at a rest-frame wavelength of 1700 A, correspond- 
ing to a star formation rate of about 15 M G yr~' (uncorrected for 
potential dust extinction). Because of the steep faint-end slopes 
of the fitted LFs, they argue that the total UV luminosity density 
(and hence star formation rate density) in this redshift range not 
dominated by X.* or even brighter gaalxies but by faint galaxies, 
mostly around luminosities of 0.1 £,*. The faint end-slopes of 
the rest-frame UV LFs beco mes even steeper at higher redshifts 
(IBouwens et al ] |2007l 120091) . Therefore, the marginal contribu- 
tion of LBGs to the total cosmic SFRD is not surprising since 
the LBG samples used in this work are based on a spectroscopic 
selection. Hence the sources must be bright enough for spec- 
troscopy, which in this cas e (e.g. for the VIMO S spectroscopic 
campaign in the CDF-S by Balestra et al. 2010) means that they 
have to have apparent magnitudes brighter than 24.5 in the B- 
and R -bands, respectively. At z = 2.0, this corresponds to an ab- 
solute B-band (so rest-frame ~ 1500 A) ma gnitude of -2 1 .5, half 
a mag nitude brighter than M* derived by ISawicki & T hompson 
(2006). At z — 3.0, the situation becomes even worse since the 
spectroscopic sample in this case is limited to objects 1.5 mag- 



nitudes brighter than M„. We therefore point out that our results 
are in very good agreement with the previous statements that 
the bulk of star formation activity takes place not in bright but 
in faint galaxies. A n independent study recently conducted by 
iTanvir et al. utilizing gamma-ray bursts (GRBs) as star 

formation rate tracers, so following an entirely different premise, 
comes to the same result. With their data for even higher red- 
shifts (z ~ 5), they argue that the bulk of star formation activity 
is not even accessible in currently available ultra-deep data sets 
such as the Hubble Ultra-Deep Fields which implies mean star 
formation rates per "typical" galaxy of less than 0.2 M Q yr~' . 



6. Summary and Conclusions 

Utilizing stacking techniques together with the newly acquired 
Chandra 4 Ms mosaic, we have investigated the X-ray luminos- 
ity and star formation activity of Lyman Break Galaxies (LBGs) 
across cosmic time. Our stacking input sample spans a redshift 
range 0.5 < z < 4.5, making use of the superb spectroscopic data 
in the Chandra Deep Field South region. Spectroscopic selection 
of LBGs guarantees a highly homogeneous sample with minimal 
contamination by interloping objects, in particular AGN contam- 
inating the measured X-ray fluxes. Our newly developed stack- 
ing algorithm, optimized for highest sensitivities in exchange for 
the loss of morphological information allows us to probe the 
LBG population down to formerly unprecedented flux density 
levels of 10- 18 mWnr 2 in the soft (0.5-2 keV) band. Our find- 
ings are summarized below: 

1. We reliably (Poissonian confidence level > 95%) detect X- 
ray emission from our LBG input sample out to redshifts of 
about 4 with a robust upper limit for z = 4.5. 

2. From the observed soft-band fluxes, we derive rest-frame 2- 
lOkeV luminosities using a photon index r = 2.0 and cor- 
recting for Galactic hydrogenabsorption. They show a nearly 
constant value with redshift, underlining the robustness of 
our stacking procedure. 

3. From the rest 2-10keV luminosities we calculated mean star 
formation rates for each redshift bin. Our results show a dis- 
tinct peak at z pe ak ~ 3.5, in good agreement with the general 
trend of various other estimates of SFR for LBGs and other 
types of galaxies. However, we point out that the various star 
formation rate indicators across the electromagnetic spec- 
trum deliver significantly different SFR values, therefore a 
thorough comparison and cross-calibration is highly needed 
also in advance of upcoming all-sky survey projects. 

4. With ancillary /T-band infrared data, we calculated stel- 
lar masses and hence specific star formation rates 
(sSFR=SFR/M„). The peaking behavior of our LBG sample 
is also seen in sSFR, underlining the widely adopted picture 
of stellar mass growth having a peak somewhere between 
redshift 2 and 4. 

5. Considering the contribution of LBGs to the total cosmic star 
formation rate density (SFRD), we find that LBGs only make 
up a tiny fraction of the total star formation activity at all in- 
vestigated redshifts, supporting the emergent notion that the 
bulk of star formation in the universe takes place in very low- 
mass, low-SFR galaxies currently escaping detection with 
all available facilities. 

All together, Lyman Break Galaxies seem to be good examples 
of "typical" galaxies at their respective redshifts but are not the 
place in the universe where most the star formation activity takes 
place. Despite their high typical star formation rates of a few to 
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several lOM yr _I at all redshifts, LBGs are just not abundant 
enough to make up a significant fraction of the total star forma- 
tion in the universe. This highly supports the proposition that the 
bulk of high-redshift star formation is going on in faint, yet un- 
detected galaxies well below X* at all redshifts z > 2 (see e.g. 
iTanvir etal.ll20ll and references therein). 

Acknowledgements. We thank our anonymous referee for his/her constructive 
comments, in particular regarding the many technical finesses of X-ray observa- 
tions. 
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